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Resumen 


Reyes-Olvera, A. L., £: Gutiérrez-González, E. (marzo-abril, 
2016). Modelación del riesgo de inundaciones en el estado 
de Tabasco en el periodo 1961-2007. Tecnología y Ciencias del 
Agua, 7(2), 99-114. 


En este trabajo se propone un indicador del riesgo de 
inundación en la región hidrográfica 30 y parte de la 29 que 
se encuentran en el estado de Tabasco, México, considerando 
tres factores: vulnerabilidad, costo y peligro. En el caso de la 
vulnerabilidad, se medirá la susceptibilidad de una región 
frente a una inundación, para esto se consideran elementos 
sociales y naturales, en donde el primero mide la disposición 
de los servicios con los que cuenta la población y el segundo 
incluye los recursos hidrológicos que atraviesan la región 
de estudio. Para el costo se considera cuánto y cómo 
afecta el daño provocado por una inundación, tomando 
en cuenta dos factores: número de habitantes y bienes que 
posee la población. En el tercer factor se consideran las 
precipitaciones máximas anuales de 1961 a 2007 tomadas de 
la base de datos Maya v. 1.0, para calcular la probabilidad de 
ocurrencia de inundación y sus periodos de retorno; para el 
peligro se propone una prueba de bondad de ajuste bootstrap, 
basada en el coeficiente de correlación muestral, estimando 
el tamaño y la potencia de la prueba usando simulación. 
Los parámetros de la distribución fueron estimados por 
máxima verosimilitud y con la prueba se verificó que las 
observaciones cumplieran con la distribución Gumbel o de 
valores extremos generalizada, elegidas por los criterios log- 
verosimilitud, AIC y correlación. 


Palabras clave: índice de riesgo, vulnerabilidad, costo y 
peligro, valores extremos, prueba bootstrap, mapa de riesgo, 
niveles y periodos de retorno, componentes principales. 


Introducción 


El agua es uno de los recursos naturales más 
valiosos de cualquier país debido a los benefi- 
cios sociales y económicos que se derivan de su 
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consciente explotación; sin embargo, junto con 
las ventajas existen también situaciones extremas, 
como inundaciones y sequías (Zucarelli, 2013). 
En el mundo, las inundaciones están aumen- 
tando más rápidamente que ningún otro desastre 
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y han causado severos problemas tanto sociales 
como económicos. Por esta razón se han desa- 
rrollado diversas investigaciones para analizar 
eventos extremos y la forma de evaluar el riesgo 
de inundación. En México, Uribe-Alcántara, 
Montes-León y García-Celis (2010) propusieron 
un mapa de índices de inundación para identifi- 
car las zonas propensas a inundaciones. 

El estudio de las inundaciones es complejo, 
pues, por ejemplo, el daño económico directo 
a un edificio inundado no sólo depende de la 
profundidad del agua sino también del uso 
del inmueble. Otros factores importantes son 
las características del edificio, las variables 
socioeconómicas y la calidad de la respuesta 
de emergencia, entre otras (Smith, 2001). En la 
Ciudad de México, Baró-Suárez, Díaz-Delgado, 
Calderón-Aragón, Esteller-Alberich y Cadena 
Vargas (2011) realizaron una investigación sobre 
el costo más probable de daños en zonas ha- 
bitacionales y propusieron una caracterización 
de costo con base en el índice de marginación 
urbana. 

Existen diferentes formas de examinar 
las inundaciones. Hunt las analiza como un 
desastre natural en el que se deben identificar 
los peligros y la vulnerabilidad para poder 
superarlos (Hunt, 2002). El Centro Nacional de 
Prevención al Desastre (Cenapred, 2006), por el 
contrario, no maneja los desastres como natu- 
rales sino como producto de las condiciones de 
vulnerabilidad y exposición. Entonces, si lo que 
se desea es prevenir los desastres, es necesario 
hablar de riesgo. Es posible señalar varias defi- 
niciones de riesgo, que se basan principalmente 
en los factores de los que depende su estima- 
ción. Además es importante mencionar que los 
efectos de las inundaciones en la sociedad no 
pueden ser representados en una cotización 
monetaria (Green, Van Der Veen, Wiertra, éz 
Penning-Rowsell, 1994). 

En algunos otros trabajos, Herrera-Díaz, 
Rodríguez-Cuevas, Couder-Castañeda y 
Gasca-Tirado (2015), así como Flowers-Cano, 
Jeffrey-Flowers y Rivera-Trejo (2014), al igual 
que Young (2002), combinan los métodos 
estadísticos y determinísticos para analizar 
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datos pasados y predecir los niveles de inun- 
dación y de ríos. Por su parte, Bayliss y Reed 
(2001) discuten diferentes formas de estudiar 
las inundaciones históricas y revisan métodos 
para incorporar datos históricos al análisis de 
frecuencias de inundaciones. Se puede incorpo- 
rar información meteorológica para fortalecer 
la estimación del riesgo; la lluvia puede ser 
un notorio indicador fiable de la inundación 
(Duncan, 2002). Sin embargo el riesgo abarca 
dos aspectos: peligro (P) y vulnerabilidad (V) 
(Mileti, 1999), aunque si bien no es la única de- 
finición que se maneja para el riesgo, sí es una 
de las más utilizadas. Ordaz incluye un tercer 
factor para la expresión del cálculo de riesgo 
por algún fenómeno natural (Ordaz, 1996), el 
costo (C) o valor de los bienes expuestos; en- 
tonces, el riesgo es producto de tres factores: 
R=VxC xP. Nótese que en realidad se trata 
de la misma definición, ya que la vulnerabilidad 
está compuesta de dos elementos: exposición 
o daño potencial, y susceptibilidad o pérdida 
(Merz, Thieken, € Gocht, 2007). Por lo tanto, 
para obtener una medida del riesgo, se debe 
trabajar sobre estos tres últimos factores. 

Debido al impacto de las inundaciones 
sobre la sociedad y las pérdidas monetarias, 
es importante la construcción de un índice de 
vulnerabilidad (Tapsell, Penning-Rowsell, € 
Tunstall, 2002). La vulnerabilidad tiene dife- 
rentes conceptos, por ejemplo, Blaikie, Wisner, 
Cannon y Davis (1994) analizan las condiciones 
socioeconómicas que aumentan el grado de 
vulnerabilidad. 

En cuanto al peligro, se deben calcular las 
probabilidades de que ocurra un fenómeno que 
pueda dañar lo expuesto (Ordaz, 1996). Por el 
efecto de las inundaciones, el riesgo puede ser 
generado por lluvias, ciclones tropicales o falla 
de obras hidráulicas. 

En este artículo, las inundaciones se estu- 
dian cuando se registran valores extremos de 
precipitaciones de lluvias y es aquí donde se 
introducirá la teoría sobre valores extremos. 

Fisher y Tippett (1928) derivaron la forma lí- 
mite de la distribución del valor máximo de una 
muestra aleatoria, encontrando las tres posibles 
distribuciones: Gumbel, Weibull y Fréchet. A 
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partir de esto se desarrollan aportaciones de in- 
terés como las de Gumbel (Gumbel, Les moments 
des distribution límite du terme maximum dune 
serie aleatoire, 1934), quien además presenta el 
primer libro de importancia (Gumbel, Statistical 
of Extremes, 1958), en el que trabajó con eventos 
extremos contemplando varias aplicaciones. Por 
su parte, Gnedenko presenta en forma general y 
con una prueba rigurosa el teorema de los tipos 
de distribución para valores extremos propuesto 
por Fisher y Tippett (Gnedenko, 1943). Otra 
aportación importante la dieron Von Mises 
(1936) y Jenkinson (1955), quienes trabajaron 
sobre la distribución de valores extremos y pro- 
pusieron una distribución que combina las tres 
familias de distribución de valores extremos, la 
cual es conocida como DVEG. De 1990 a 2011, 
el desarrollo estadístico ha sido sustancial en 
el área, como en el caso, entre otros, de Tawn 
(1992), Rosbjerg y Madsen (1996), Coles y Dixon 
(1999, Likelihood-Based Inference for Extreme Values 
Models Extremes) y Coles (2001, An Introduction 
to Statistical Modeling of Extreme Values). 

Obtenida una medida del riesgo, se im- 
plementan los mapas de peligro y riesgo por 
inundación para ubicar las zonas en diferentes 
categorías de riesgo y desarrollar un plan de 
manejo. Tawatchi y Mohammed (2005), y Leh- 
ner y Dóll (2001) proporcionan un ejemplo para 
aproximar a escala global el análisis con mapas 
concernientes a la situación de inundación en 
Europa; en cuanto a México, se tienen los mapas 
de riesgo a escala municipal por inundaciones 
(Cenapred, 2006). 

En México, en la región del sistema Grijalva- 
Usumacinta del estado de Tabasco se han 
presentado inundaciones en los últimos años. 
Una de ellas, en 2007, fue una de las de mayor 
intensidad en los últimos 50 años (Rivera-Trejo, 
Soto-Cortés, €: Barajas-Fernández, 2009). Por su 
parte, Arreguín-Cortés, Rubio-Gutiérrez, Do- 
mínguez-Mora y Luna-Cruz (2014) analizaron 
los factores que influyen en las inundaciones 
de la planicie tabasqueña, como la ausencia de 
ordenamiento territorial adecuado y la defores- 
tación de la parte alta de las cuencas. 


Por esta razón, en la presente investigación 
se calculan los índices para el riesgo de inunda- 
ción y sus periodos de retorno con cuatro esce- 
narios de peligro: 100, 150, 200 y 250 mm. Para 
esto se propone una prueba de bondad de ajuste 
bootstrap paramétrica basada en el coeficiente de 
correlación muestral, y se comprueba el tamaño 
y potencia de la prueba por simulación realizan- 
do 10 000 réplicas. 


Objetivo general 


Modelar el riesgo de inundaciones en el estado 
de Tabasco por causas climáticas extremas. 


Objetivos específicos 


+ Proponer un modelo estadístico basado en 
la teoría de valores extremos que permita 
determinar el comportamiento y las 
tendencias de las precipitaciones extremas. 

e  Calibrar el modelo estadístico usando una 
base de datos histórica de precipitación y 
niveles hidrométricos del estado de Tabasco. 

e Construir un mapa de riesgo de 
inundaciones por municipio en el estado 
de Tabasco. 


Metodología 


Por lo común se cree que el peligro es el único 
responsable de los desastres, pero en buena 
medida es la sociedad la que se expone con su 
infraestructura, organización y cultura para 
enfrentar fenómenos con alto grado de peligro. 
Por lo tanto, se debe entender que los desastres 
no sólo son naturales, sino que también son 
producto de condiciones de vulnerabilidad y 
exposición (Cenapred, 2006). 

Por lo anterior, es necesario encontrar una 
medida que represente el riesgo y sus cambios 
a través del tiempo; puede ser estadísticamente 
difícil de modelar debido a la gran cantidad de 
variables que intervienen en el modelo. Una 
forma de tratar y aligerar la solución a este 
problema puede ser introduciendo los núme- 
ros índice. Los índices, desde su aparición, han 
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llegado a ser cada vez más importantes como 
indicadores de los diferentes cambios que ro- 
dean a la sociedad. 

En este trabajo se van a construir los índices 
y mapas de riesgo a escala municipal en el 
estado de Tabasco. Para el riesgo se utilizará 
la definición de Ordaz (1996), quien maneja al 
riesgo como el producto de tres factores (ver ex- 
presión (1)). Para cada uno de éstos se calculará 
un indicador: 


R=VxCxP (1) 


En donde V representa la vulnerabilidad, C, 
el valor de los bienes expuestos y P es el peligro 
o probabilidad de ocurrencia de un fenómeno 
que pueda dañar lo expuesto. 


Vulnerabilidad 


La vulnerabilidad puede tener diferentes 
conceptos. En este artículo se considera como 
una medida de qué tan susceptible es un bien 
expuesto a la ocurrencia de un fenómeno 
perturbador (Ordaz, 1996). Para su análisis se 
consideran dos variables: por un lado, los ser- 
vicios con los que cuenta la comunidad (energía 
eléctrica, agua potable y drenaje) y, por el otro, 
sus principales recursos hidrológicos (ríos y 
lagunas). 

En el caso de los servicios, se construye un 
índice de servicios a través de tres variables: 


PO_NDEE: porcentaje de ocupantes que no 
dispone de energía eléctrica. 

PO_NDAE: porcentaje de ocupantes que no 
dispone de agua entubada. 

PO_NDD: porcentaje de ocupantes que no 
dispone de drenaje. 


Estos factores son ponderados mediante el 
uso de componentes principales (CP) del aná- 
lisis multivariado. El uso de CP garantiza una 
combinación lineal con máxima varianza que 
brinda los mejores coeficientes y que explica de 
modo adecuado la situación de la comunidad. 

La segunda variable, recursos hidrológicos 
con los que cuenta cada municipio, es considera- 
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da debido a que la región de Tabasco está a muy 
poca altura sobre el nivel del mar; es una zona 
de ríos, lagunas y pantanos. Además, a lo largo 
de su historia, Tabasco ha sufrido inundaciones 
por el desbordamiento de los ríos que atraviesan 
la entidad. Por ejemplo, la inundación de 2007 
(Rivera-Trejo et al., 2009) en el municipio de Pa- 
raíso se debió al exceso de lluvia que coincidió 
con el oleaje en el Golfo de México, que evitó las 
descargas de los ríos Grijalva y Carrizal hacia el 
mar. Entonces, para los recursos hidrológicos, 
se construye un índice de ríos a través de una 
ponderación por municipio que depende del 
número de ríos que se encuentren en él, además 
de considerar la capacidad de cada uno de los 
ríos según la información del nivel de agua 
máximo ordinario (NAMO). 


Costo de los bienes expuestos 


La componente costo, C, de la ecuación (1) 
mide la cuantía de lo que es susceptible de 
afectarse durante la ocurrencia de un fenóme- 
no perturbador (Ordaz, 1996). Entonces, para 
la construcción de un indicador de costo, se 
consideran los bienes con los que cuenta la 
población (televisión, lavadora, refrigerador, 
computadora, automóvil) y el número de habi- 
tantes por comunidad. El número de habitantes 
se toma en cuenta por el supuesto de que “A 
mayor población, mayor riqueza”. Entonces, el 
índice de costo de bienes expuestos se construye 
con base en los siguientes cinco indicadores de 
bienes, que después serán ponderados mediante 
el uso de componentes principales: 


PODT: porcentaje de ocupantes que dispone de 
televisión. 

PODR: porcentaje de ocupantes que dispone de 
refrigerador. 

PODL: porcentaje de ocupantes que dispone de 
lavadora. 

PODC: porcentaje de ocupantes que dispone de 
computadora. 

PODA: porcentaje de ocupantes que dispone de 
automóvil o camioneta. 
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Peligro 


El peligro, P, de inundación en el contexto de es- 
ta investigación se define como la probabilidad 
de ocurrencia de una situación de inundación 
potencialmente dañina en un área dada y con 
cierto periodo especificado. Para el análisis 
de este componente, la información a utilizar 
dependerá del tipo de inundación que se analice 
y pueden definirse de acuerdo con: 


e La duración de la inundación: lenta o súbita. 

e El mecanismo que las genere: pluviales, 
fluviales, costeras o falla de obras 
hidráulicas. 


En este caso, se calcula un valor para el 
peligro tomando a la lluvia como un indicador 
fiable de la inundación; no se estudian las 
inundaciones dadas directamente por fallas de 
obras hidráulicas, aunque es importante hacer 
notar que cuando se presentan valores extremos 
de precipitación es muy probable que ocurran 
fallas en obras hidráulicas. Por tanto, para el cál- 
culo del factor peligro se analizan los datos de 
precipitación en Tabasco y zonas cercanas; estas 
últimas se consideran, pues las inundaciones no 
dependen sólo de la precipitación local, sino de 
la precipitación en la cuenca aguas arriba del 
punto de desborde. 


Datos de precipitación 


La información de la precipitación se obtuvo 
de las mediciones diarias —desde 1961 hasta 
2007— de precipitación tomadas en las estacio- 
nes registradas por el Servicio Meteorológico 
Nacional (SMN) con la herramienta CLICOM, 
en su versión interpolada en malla regular de 
0.2 grados de longitud por 0.2 grados de latitud 
(Maya v. 1.0). 

Para obtener la base de datos de cada mu- 
nicipio se delimitó una región para cada uno. 
Esto se realizó con ayuda de la información de 
las subcuencas, cuyo caudal influye en el mu- 
nicipio de estudio (ya sea parcial o totalmente). 
En caso de que parte de alguna subcuenca entre 


en el municipio, pero la aportación de su cau- 
dal no influya, se toman sólo los nodos que se 
encuentren dentro del municipio. Con esto se 
determina cuáles son los nodos que influyen en 
el análisis, así se tendrán tantas series de datos 
como el número de nodos incorporados en el 
estudio. Entonces, para obtener una sola serie, 
se realiza una combinación lineal, que resulte 
adecuada, de los nodos considerados. Para esto 
nuevamente se utiliza la teoría de componen- 
tes principales; en este caso, las variables son 
los nodos del municipio. Las CP obtenidas 
corresponden a los pesos dados a cada nodo 
para obtener la combinación lineal deseada que 
finalmente será la serie de datos que represente 
al municipio. 


Análisis de los datos 


Dadas las series de datos de cada municipio, 
se extraen las precipitaciones máximas anua- 
les. Sea X,, Xp». X, Una muestra de la variable 
aleatoria X:: precipitación en el municipio 1. 
Entonces se forman los vectores x= (x1, ..., x!), 
x= (12, ..., 15), 1% = (xt, ..., 1). Con los vec- 
tores se calculan los máximos anuales, 
Ma, = máxlx;, ..., Xx), 1 =1, ..., k, que gracias a 
la teoría de valores extremos se consideran ¿id 
(independientes e idénticamente distribuidos), 
entonces F(x) = P(M < x) = [E(1)Y. Por el teo- 
rema de Fisher y Tippett (1928) se asume que 
M, sigue una de las distribuciones de valores 
extremos (DVE) mostradas en (2), (3) y (4): 


Gumbol F(2)xp exp (Es) ! xEeR (2) 
o 
0 x<u 
Fréchet: F(x) = Pr (3) 
exp (5) x>u 
o 
0 x<u 
Weibull: F(x) = Ev (4) 
1-exp E) x>u 
o 
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con parámetros 4 ER, 0 >0 y e >0 o bien consu 
generalización dada por Gnedenko, la distribu- 
ción de valores extremos generalizada (DVEG), 


a] (5) 
o 5 


donde y, = máx(y, 0) JER,o>0yeER. 
En esta investigación se estudió el ajuste con 


mostrada en (5): 


la distribución Gumbel, por ser una de las más 
empleadas en el campo de las inundaciones en 
México, y la DVEG, cuyo uso se ha ampliado en 
los últimos años, pero aún no es muy aplicada 
para las precipitaciones máximas en México. 
Los parámetros u, O y e se estiman por el méto- 
do de máxima verosimilitud, para aprovechar 
las propiedades asintóticas de estos estimadores 
y Obtener los intervalos de confianza de forma 
sencilla. La función log-verosimilitud para la 
DVEG y la distribución Gumbel están dadas en 


(6): 
te to 


-.tog(0)-(2+1))os 


i=l 


1/8 


sig0 


l(u,0,e;x) = Nte > 


-nlog(0) - Sexp 
i=1 


=x.- ul al-x.-u|. 
—— 1 -+| 2 1 sie=0 
AA 


(6) 


Sean (ñ, 0,£) los estimadores de máxima ve- 
rosimilitud (EMV) de (u,o,e), éstos se calculan 
con las librerías evd y VGAM del Proyecto R. 
Para la estimación por intervalos, se usa el 
hecho de que los EMV se distribuyen asintóti- 
camente normal ÓN, (0,5 (6) donde 1(8) es 
la matriz de información de Fisher. Entonces, el 
intervalo de confianza a (1 — a1)100% para cada 
componente de 6 está dado por 0,Zo) Cas 
en donde Z,,, es el valor de la variable nor- 
mal estándar con área derecha igual a 0/2, 
y c, es el ¡-ésimo elemento de la diagonal 17 (6). 
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Elección del modelo 


La elección del modelo se basó principalmente 
en el criterio de información de Akaike (AIC), 
AIC(k) = 2log(L(0,)) + 2k, donde k es el núme- 
ro de parámetros del modelo y log(L(0,)) es el 
logaritmo de la verosimilitud del modelo, con- 
siderando mejor modelo si tiene menor valor 
de AIC. Un segundo criterio para modelos con 
la misma cantidad de parámetros es el valor 
del logaritmo de la verosimilitud (log-veros), 
eligiendo el modelo que tenga mayor valor de 
log-verosimilitud. El tercer criterio asume que 
existe una relación lineal entre la distribución 
empírica y la distribución teórica correspon- 
diente, Gumbel (2) o VEG (5), por lo que se 
emplea el coeficiente de correlación muestral 
(7) como medida de asociación y se prefiere el 
modelo cuyo coeficiente de correlación muestral 
sea más cercano a uno: 


Prueba bootstrap basada en el coeficiente 
de correlación muestral 


Sea X,..., X, una muestra aleatoria de una 
función de distribución F(x) = P (X < x). El 
contraste de hipótesis a probar está dado por: 


Hi EF e) 0. H:fle EF" le) 18) 


donde F*(x) denota una familia de densidades 
específica; en este caso se trata de la función de 
densidad Gumbel y la DVEG. 

Las densidades que se analizan en esta 
investigación son del tipo localidad y escala 
(Gutiérrez-González, Panteleeva, e Córdoba- 
Lobo, 2012) cuando se fija un valor para el pa- 
rámetro de forma. Así, con base en la invarianza 
bajo la transformación de los parámetros de 
localidad y escala, se propone una prueba de 
bondad de ajuste bootstrap paramétrico basada 
en el coeficiente de correlación muestral, como 
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una estadística de prueba. Una prueba similar 
a ésta para la distribución log-gamma genera- 
lizada fue presentada por Gutiérrez-González, 
Villaseñor-Alva, Panteleeva y Vaquera-Huerta 
(2013). 

Si H, de (8) se cumple, entonces la distribu- 


ción de f(x) es tal que F(x;u,0,8) nt). 
o 
Entonces, para una variable aleatoria X con 


función de densidad f(x), se tiene que la función 
F (x;2) le corresponde su función inversa F-Xy;e) 
eny. 

Sean Xy X) las estadísticas de orden de 
una variable aleatoria, entonces la función de 
distribución empírica correspondiente está dada 
por: 


+05. 
El )= > Si X(] SEX <Xf5ay 
1 
l si x> Xi) 
donde X¿, es el valor de la ¡-ésima estadística 
de orden. 


Por el teorema de Glivenko-Cantelli se es- 
tablece una convergencia casi segura, cuando 
n => oo, entre la distribución empírica F, (x) y la 
distribución teórica, que puede ser F(x) o G(x); 
por lo tanto, es factible establecer para la DVEG 
una aproximación: 


G,(x)=G(x)=> : [-log[G, (1 ay ll 


Similarmente, para la distribución Gumbel 
resulta: 


09 =rt) ts relata) E) 


Ambas expresiones se pueden reescribir 
como Y = f + aZ y se espera que bajo H, exista 
una fuerte relación lineal entre Y y Z cuando se 
sustituye al parámetro e por un estimador con- 
sistente (2). Entonces, para probar el contraste 
de hipótesis de (8), se propone como estadística 
de prueba al coeficiente de correlación muestral 


dado en (7), en el cual se espera que bajo H, 
la distribución de r,(Y,Z) estará concentrada 
en las proximidades de 1. Por tanto, la regla de 
decisión para una prueba de tamaño a, con a € 
(0,1), un valor conocido, consiste en rechazar H, 
cuando r,(Y,Z) <r... 

La prueba bootstrap con base en los EMV se 
realiza siguiendo los siguientes pasos: 


1. Dadas las observaciones x.,... x,, calcular un 
valor para el estimador de e, denotado por 
Eb: 

2. Con las observaciones se calcula su coefi- 
ciente de correlación muestral 7. 

3. A partir de £, iniciar un ciclo bootstrap. 

a) Generar una muestra bootstrap de ta- 
maño n de la DVEG o la distribución 
Gumbel, según el caso. 

b) Con la muestra del inciso (3a) calcular 
una estimación para el parámetro de 
forma, denotándola por £,. 

c) Cong, del inciso (3b) se calcula el coefi- 
ciente de correlación muestral 7. 

4. El ciclo de bootstrap del inciso 3 se repite 
m veces para calcular, 7,,..., ?,, en donde m- 
cantidad de estimaciones de bootstrap, para 
determinar el cuantil bootstrap. 

5. Con el ciclo terminado, ?, ..., F, son ordena- 
dos en forma no decreciente, denotándolos 
con 7. Entonces, Í, sf, <---<?,, y Obtener el 
cuantil a,, sea éste F.. 

6. Regla de decisión. Comparar r,, con el cuan- 
til a del paso 5. 

a) Sir, <?, se rechaza H, al nivel de signi- 
ficancia a. 

b) Si r, >f,, no se rechaza H, al nivel de 
significancia a. 


Antes de iniciar con las pruebas, se determi- 
na si éstas conservan el tamaño nominal de la 
prueba y presentan buenas potencias. 

Los resultados de los tamaños de la prueba 
para 10 000 réplicas de ambas distribuciones se 
muestran en el cuadro 1. Se puede concluir que 
la prueba bootstrap conserva el tamaño nominal 
de la prueba para las distribuciones Gumbel y 
VEG. 
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Cuadro 1. Tamaños de prueba para la DVEG y la distribución Gumbel. 


DVEG con e =0.1 DVEG con e =0.5 Distribución Gumbel 


n | 0.010 | 0.025 | 0.050 | 0.075 | 0.100 | 0.010 | 0.025 | 0.050 | 0.075 | 0.100 | 0.010 | 0.025 | 0.050 | 0.075 | 0.100 


30 | 0.012 | 0.026 | 0.040 | 0.060 | 0.086 | 0.004 | 0.011 | 0.039 | 0.059 | 0.084 | 0.008 | 0.023 | 0.039 | 0.060 | 0.082 


50 | 0.008 | 0.016 | 0.030 | 0.046 | 0.064 | 0.004 | 0.012 | 0.038 | 0.063 | 0.095 | 0.010 | 0.030 | 0.055 | 0.080 | 0.100 


75 | 0.002 | 0.010 | 0.020 | 0.036 | 0.050 | 0.006 | 0.018 | 0.038 | 0.058 | 0.086 | 0.012 | 0.021 | 0.043 | 0.069 | 0.097 


100 | 0.004 | 0.016 | 0.026 | 0.046 | 0.062 | 0.006 | 0.018 | 0.041 | 0.059 | 0.085 | 0.010 | 0.018 | 0.046 | 0.078 | 0.100 


150 | 0.004 | 0.010 | 0.044 | 0.058 | 0.070 | 0.007 | 0.020 | 0.046 | 0.069 | 0.096 | 0.011 | 0.029 | 0.053 | 0.077 | 0.104 


En el caso de la potencia de la prueba se uti- 
lizan tres diferentes distribuciones alternativas: 
normal estándar, t-Student con 4 gl y Cauchy 
(0.2). Los resultados para 10 000 réplicas se 
muestran en el cuadro 2. 


Periodos de retorno 


Para el análisis del factor peligro es de interés 
conocer el tiempo promedio de espera para que 


ocurra una lluvia extrema que pudiera causar 
algún problema de inundación. El periodo de 
retorno es el tiempo medio para que se repita 
un evento semejante. Entonces sea X,, X,, ... 
una sucesión de variables aleatorias ¡id con 
función de distribución continua F y u ER. 
Considerando (X, > uj un evento extremo, se 
define el tiempo de la primera excedencia como 
T(u) = mínfi > 1:X, > u). Con esto, el tiempo 
promedio de espera para que ocurra un evento 


Cuadro 2. Potencia de la prueba para la DVEG. 


Alternativa 
' Ñ Normal (0.1) t-Student (4) Cauchy (0.2) 

0.010 0.376 0.398 0.658 

A 0.050 0.662 0.610 0.824 
0.075 0.734 0.660 0.862 

0.100 0.774 0.704 0.882 

0.010 0.534 0.396 0.884 

0.050 0.850 0.704 0.964 

ss 0.075 0.900 0.782 0.976 
0.100 0.924 0.814 0.981 

0.010 0.536 0.420 0.960 

0.050 0.930 0.790 0.996 

Ús 0.075 0.956 0.856 0.998 
0.100 0.964 0.892 1.000 

0.010 0.552 0.604 1.000 

00 0.050 0.976 0.862 1.000 
0.075 0.992 0.894 1.000 

0.100 0.992 0.930 1.000 

0.010 0.574 0.781 1.000 

0.050 0.994 0.966 1.000 

ed 0.075 0.998 0.987 1.000 
0.100 1.000 0.992 1.000 


+» ISSN 0187-8336 


Reyes-Olvera y Gutiérrez-González, Modelación del riesgo de inundaciones en el estado de Tabasco en el periodo 1961-2007 


extremo es E[T(u). Si la probabilidad del evento 
[(X * > u) es p, su periodo de retorno es p. En- 
tonces, usando la distribución de probabilidad 
acumulada, para el evento (X, > u), el periodo 
de retorno se calcula con (9): 


== (9) 


En específico, el nivel de retorno asociado 
con el periodo de retorno p” es el (1 — p)-ésimo 
cuantil de la distribución. En el caso de la dis- 
tribución Gumbel y la de VEG, el cuantil 1 — p 
está dado en (10): 


ro] -e : 
o u-S|1-(1og(1=p)) si 2x0 de 


u+0|-log(-log(1-p))] sig=0 


Los cuantiles se calculan a partir de las esti- 
maciones de los parámetros. 


Resultados y discusión 


El índice de riesgo se calcula con la fórmula (1). 
Para los factores de vulnerabilidad y de costo se 
construyen dos indicadores con base en la dis- 
posición de servicios, de los bienes con los que 
cuenta la población y del número de habitantes 
por municipio. La información fue tomada de 
la base de datos de INEGI, que corresponde al 
censo de población y vivienda 2010. 


Índice de vulnerabilidad 


La vulnerabilidad se compone de dos factores: 
social y natural. El factor social se refiere al indi- 
cador de servicios (Ind, ) que se obtiene a partir 
de una combinación lineal de las tres variables 
PO_NDEE, PO_NDAE y PO_NDD. Los pesos 
de la combinación lineal se obtienen por CP. Con 
la información de las tres variables se calcula 
la matriz de covarianzas y de ésta, los valores 


propios de cada una de las variables, resultando 


290.367, 6.808 y 1.228, respectivamente. Después 
se ponderan los valores propios. Resultan los 
pesos de cada una de las tres variables. Para 
obtener el Ind... , los resultados se muestran en 
el cuadro 3. 

Para el factor natural se propone el indicador 
ríos (In), en donde se consideran los princi- 
pales recursos hidrológicos con los que cuenta 
cada municipio. El indicador se construye pon- 
derando el número de ríos que forman parte del 
municipio y la capacidad que tiene cada uno, 
medida por el nivel de agua máximo ordinario. 
Los resultados del In. se muestran en el cuadro 
3. 

Por último, para obtener el indicador de 
vulnerabilidad (In, ), se ponderaron ambos 
indicadores, In... e In,,.. Los resultados se 
muestran en el cuadro 3. 


Índice de costos 


Para el análisis del factor de costo intervienen 
dos variables de interés. En primer lugar, se 
tiene el logaritmo del número de habitantes en 
cada municipio, esto basado principalmente en 
el supuesto de que a mayor población mayor 
riqueza. Se toma el logaritmo de la población 
(log,.») sólo para mantener una escala más 
apropiada. Los resultados se muestran en el 
cuadro 4. 

La segunda variable que se toma en cuenta 
es un indicador de los bienes (In,,....) con los que 
cuenta la población del municipio en estudio. 
Después, en la construcción de este indicador se 
consideraron cinco bienes que estarían expues- 
tos a sufrir cierto daño ante una inundación y 
se denotan de la siguiente manera: televisión 
(PODT), refrigerador (PODR), lavadora (PODL), 
computadora (PODC), y automóvil o camioneta 
(PODA). En el cálculo del indicador de bienes 
se consideraron los porcentajes de la población 
que cuenta con tal o cual bien. 

Para obtener una combinación lineal adecua- 
da de estas cinco variables se realiza un análisis 
de CP para generar los pesos adecuados de 
cada variable. De forma similar al caso de los 
servicios, se calcula la matriz de covarianzas y 
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Cuadro 3. Índice de vulnerabilidad. 


Núm. Municipios In... In... Mon Categoría 
1 Balancán 0.282 0.128 0.221 Medio 
2 Cárdenas 0.403 0.785 0.556 Muy alto 
3 Centla 1.000 0.161 0.664 Muy alto 
4 Centro 0.000 0.957 0.526 Muy alto 
5 Comalcalco 0.242 0.552 0.366 Alto 
6 Cunduacán 0.531 0.785 0.632 Muy alto 
7 E. Zapata 0.001 0.161 0.065 Bajo 
8 Huimanguillo 0.810 0.101 0.526 Muy alto 
9 Jalapa 0.123 0.827 0.405 Alto 
10 Jalpa 0.166 0.684 0.373 Alto 
11 Jonuta 0.123 0.161 0.138 Bajo 
12 Macuspana 0.496 0.665 0.564 Muy alto 
13 Nacajuca 0.016 0.392 0.166 Bajo 
14 Paraíso 0.112 1.000 0.467 Alto 
15 Tacotalpa 0.367 0.000 0.220 Medio 
16 Teapa 0.166 0.060 0.124 Bajo 
17 Tenosique 0.256 0.161 0.218 Bajo 
Cuadro 4. Índice de costo. 
Núm Municipios 108. Mas Mo. 
1 Balancán 4.754 0.274 2.514 
2 Cárdenas 5590 0.514 2.955 
= 3 Centla 5.009 0.238 2.624 
3 4 Centro 5.806 1.000 3.403 
E 5 Comalcalco 5.285 0.447 2.866 
E 6 Cunduacán 5.102 0.412 2 
3 Z E. Zapata 4.470 0.672 2.571 
E 8 Huimanguillo 5.204 0.286 2.770 
= 9 Jalapa 4.561 0.593 2.577 
al 10 Jalpa 4.921 0.546 2.734 
E 11 Jonuta 4.470 0.156 2.313 
> 12 Macuspana 5.185 0.428 2.806 
E 13 Nacajuca 5.061 0.742 2.901 
= 14 Paraíso 4.938 0.702 2.820 
E 15 Tacotalpa 4.666 0.000 2.333 
3 16 Teapa 4.729 0.400 2.564 
E 17 Tenosique 4.771 0.421 2.596 


Te 
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se calculan los valores propios para las cinco va- 
riables mencionadas arriba, obteniendo 243.062, 
17.619, 5.178, 3.471 y 1.405, respectivamente. 
Realizando una ponderación a partir de los va- 
lores propios, resultan los pesos de cada una de 
las cinco variables. Los resultados se muestran 
en el cuadro 4. 

Por último, para obtener el indicador de 
costo (In__..), se hace una ponderación de los 
factores log, ,, y My; que influyen en este indi- 
cador. Cabe mencionar que no existe alguna ra- 


costo 


bienes 


zón relevante para darle mayor peso a un factor 
en particular, entonces se restringe a ponderar 
con el mismo peso a ambos factores, es decir, 
0.5 para los dos. Los resultados se muestran en 
el cuadro 4. 


Factor peligro 
Obtención de la información 


Para calcular los índices de peligro, la in- 
formación de precipitación se obtuvo de las 
mediciones proporcionadas en la base de 
datos que el Servicio Meteorológico Nacional 
maneja con la herramienta CLICOM, desde 
1961 hasta 2000, en su versión interpolada en 
malla regular de 0.2 grados de longitud por 0.2 
grados de latitud (la denominada Maya v. 1.0). 
La información proviene de algo más de 5 000 
estaciones climatológicas de la base de datos. La 
información más reciente (2001-2007) proviene 
de las mediciones diarias del subconjunto de 
estaciones climatológicas que reporta en tiempo 
casi real; estos datos fueron obtenidos con base 
en su longitud y latitud. 

Los datos utilizados son los de precipitación 
pluvial y en la presentación de nodo por nodo. 
Lo que se tiene es una lista de 4 542 archivos de 
texto en donde cada uno de éstos corresponde 
a un nodo y cada archivo contiene una serie de 
datos de precipitación pluvial desde enero de 
1961 hasta diciembre de 2000. Además de esta 
información, se conoce la precipitación de 2001 a 
2007 de las estaciones existentes en la república 
mexicana. Se cuenta con un archivo por año con 
cada uno con los siguientes campos: latitud, lon- 
gitud, clave y nombre de la estación, altura, mes 


y días. Para aprovechar la información extra se 
trabajó con las estaciones que rodean al nodo. 
Con cada nodo que se encuentra involucrado 
en el análisis se calcula un valor máximo de 
precipitación de 2001 a 2007 y se obtiene la pre- 
cipitación que corresponde al nodo de interés. 

Los nodos de interés son los que intervienen 
en cada municipio. Para esto se localizan las 
cuencas que conforman al estado de Tabasco. La 
mayor parte de la superficie del estado (75.22%) 
se ubica en la Región Hidrográfica número 30, 
o Región del Sistema Grijalva-Usumacinta, for- 
mada por las cuencas hidrográficas del Grijalva, 
Usumacinta y de la laguna de Términos, que 
en Tabasco ocupan el 41.45, 29.24 y 4.53%, res- 
pectivamente. Mientras que el 24.78% restante 
se encuentra dentro de la Región Hidrográfica 
número 29 o Región del Coatzacoalcos, formada 
por dos cuencas: la del Coatzacoalcos, Tonalá 
y lagunas del Carmen, y la Machona, siendo 
esta última la única que tiene representación en 
el estado. Ambas regiones están consideradas 
como las más húmedas del país, en primer y 
segundo lugar, respectivamente. 

Con base en la información que se conoce de 
vulnerabilidad y costo se trabaja por municipio. 
Entonces, para obtener los datos de precipita- 
ción en cada municipio se utilizó la información 
de las subcuencas que forman parte de cada uno 
de éstos, ya sea que se encuentre completamente 
solo o como parte de ella, y se toman en cuenta 
aquellas cuya aportación de su caudal influye 
en el municipio de estudio. 


Análisis de la información 


En el análisis de los datos se requiere de una 
serie de observaciones para cada municipio, 
pero en este caso se tienen tantas series como 
el número de nodos que se involucren. Enton- 
ces, para obtener una sola serie se realiza una 
combinación lineal que resulte adecuada de los 
nodos que entran en el municipio. Para tener la 
mejor combinación lineal, se utiliza CP con los 
nodos de cada municipio, esto proporcionará 
la serie final de los datos de precipitación del 
municipio de interés. 
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Para el factor peligro, a los datos de preci- 
pitación se les extrajeron los máximos anuales 
a los que se ajustaron la distribución Gumbel 
y de VEG; en ambas se realizó una estimación 
puntual y por intervalos. Se consideraron los 
tres criterios mencionados en la sección anterior 
y se eligió el modelo en el que dos o más crite- 
rios le favorecieran para modelar el municipio 
de interés. 

En los cuadros 5 y 6 se muestran los re- 
sultados de las estimaciones y en la figura 1, 
las gráficas correspondientes al municipio de 
Balancán, Tabasco. La misma estimación se 
realizó para los 16 municipios restantes. En el 
cuadro 7 se ejemplifican con cuatro municipios 
los modelos elegidos. 


Una vez obtenidos los modelos para los 17 
municipios, se calculan sus respectivos periodos 
de retorno. En el cuadro 8 se muestran los re- 
sultados bajo cuatro escenarios de peligro: 100, 
150, 200 y 250 mm. Por ejemplo, para Balancán 
se tiene que el tiempo promedio para que ocurra 
una lluvia mayor o igual a 150 mm es de 3.5 
años. 

Finalmente, con los resultados de los tres 
factores se construye el Índice de Riesgo de 
Inundación, en el peligro se toma la probabili- 
dad p del evento [X, > x), donde x = 150 mm, ya 
que la experiencia sugiere que a partir de este 
valor se han presentado problemas de inunda- 
ción. Los resultados se muestran en el cuadro 
9 y en la figura 2 se tiene el mapa de riesgo de 
inundación. 


Cuadro 5. Estimación de parámetros. 


Estimación de parámetros, Balancán 


Distribución Escala Forma Localidad Log-veros AIC Correlación 
VEG 27.840 0.129 117.294 -234.053 474.106 0.974 
Gumbel 29.368 119.298 -234.831 473.662 0.972 
Cuadro 6. Intervalos de confianza. 
IC VEG Gumbel 
Escala (21.002, 34.678) (27.172, 31.564) 
Forma (-0.092, 0.351) 
Localidad (108.315, 126.274) (110.456, 128.139) 
Cuadro 7. Modelos elegidos por municipio. 
Municipio Modelo elegido 
a 
x-117.2941| 012 
á F(x)=exp3-|1+0.129| == 
Balancán ( ) p 27.84 ] 
a G(x) mej - 142.057 
ra 30.375 — 
a m G(x) =expl-exp|- [E x-163.261 
uimanguillo 753.403 
PUN 
x-91.188 | 0282 
j F(x)=expx-|1+0.282| === 
Tenosique ( ) Pp | | 29.447 ] 
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[19] 
o S Frecuencia relativa 
A Precipitación en mm S 
o 
S Sy 
S S — VEG 
Gumbel 
LO 
B S 
A S 
3 o 
Ñ 3 
o A ———— E E 
100 150 200 250 
1960 1970 1980 1990 2000 Año Precipitación en mm 
Dispersión de la precipitación en mm Ajuste de las distribuciones 
Figura 1. Gráficas de la estimación para el municipio de Balancán. 
Cuadro 8. Periodos de retorno por municipio. 
E Periodos de retorno 
Municipio 
100 mm 150 mm 200 mm 250 mm 
Balancán 1.174 3.513 12.860 41.492 
Cárdenas 1.206 3.370 14.765 72.082 
Centla 1.094 2.674 8.852 27.072 
Centro 1.019 1.862 7.249 35.443 
Comalcalco 1.437 4.150 13.180 38.315 
Cunduacán 1.241 2.954 10.099 38.515 
E. Zapata 1.038 2.077 6.397 19.303 
SS 
Huimanguillo 1.035 1.359 2.460 5.405 > 
mo 
Jalapa 1.018 1.745 6.122 27.262 pa 
Jalpa 1.109 2.328 7.870 30.891 les 
Jonuta 1.112 3.104 11.265 36.163 eS 
vu 
Macuspana 1.205 3.364 12.205 41.889 3 
Nacajuca 1.112 2.303 7.594 29.071 E 
Paraíso 1.049 2.843 16.789 115.446 3 
Tacotalpa 1.065 1.410 2.518 5.785 x 
Teapa 1.054 1.434 2.623 5.702 5 
Tenosique 1.295 2.191 13.073 27.014 =] 
? 
Discusión de resultados la hace una mejor prueba para las DVEG que 5, 
<= 
la presentada por Campos-Aranda (2001). En 3 
La prueba de bondad de ajuste boostrap propues- investigaciones previas sobre lluvias se tienen 3 
ta conserva los valores nominales del tamaño de trabajos que abarcan sólo los valores históri- 5 
5 2 4 - >= 
la prueba y las potencias bastante altas, lo que cos (González-Camacho, Pérez-Rodríguez, éz z 
S 
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Cuadro 9. Índice de riesgo propuesto. 


Núm. Municipios C VA P150 IR Categoría 
1 Balancán 2.514 0.221 0.285 0.158 Bajo 
2 Cárdenas LOS 0.556 0.290 0.487 Medio 
3 Centla 2.624 0.664 0.374 0.652 Alto 
4 Centro 3.403 0.526 0.345 0.962 Muy alto 
5 Comalcalco 2.866 0.366 0.241 0.253 Medio 
6 Cunduacán 2.757 0.632 0.316 0.590 Alto 
7 E. Zapata 2.571 0.065 0.481 0.080 Bajo 
8 Huimanguillo 2.770 0.526 0.718 1.000 Muy alto 
9 Jalapa 2.577 0.405 0.555 0.580 Alto 
10 Jalpa 2.734 0.373 0.421 0.438 Medio 
1 Jonuta 2.313 0.138 0.322 0.103 Bajo 
12 Macuspana 2.806 0.564 0.297 0.470 Medio 
13 Nacajuca 2.901 0.166 0.425 0.210 Bajo 
14 Paraíso 2.820 0.467 0.241 0.463 Medio 
15 Tacotalpa 2.333 0.220 0.682 0.364 Medio 
16 Teapa 2.564 0.124 0.680 0.221 Bajo 
17 Tenosique 2.596 0.218 0.481 0.258 Medio 


Índice de riesgo 


HN V uy ato (1.00-0.76) 
MA Arto (0.75-0.51) 
[ET] Medio (0.50-0.26) 
|] Bejo (0.25-0.00) 


Figura 2. Mapa del índice de riesgo propuesto. 


Ruelle, 2011), (Coronel-Brizio € Llanos-Arias, habitacionales (Baró-Suárez, Díaz-Delgado, 
1996; Díaz-Delgado, Bá, é Trujillo-Flores, 1999); Calderón-Aragón, Esteller-Alberich, € Cadena- 
similarmente, para la vulnerabilidad existen tra- Vargas, 2011). En la investigación hecha para el 
bajos que se relacionan con este factor (Green et riesgo de inundaciones, el índice propuesto está 
al., 1994) y lo mismo para el de costos en zonas basado en el ACP e incluye todas las variables 
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sobre las que se tiene información, agrupadas 
en los tres factores mencionados (Ordaz, 1996), 
proporcionando de esta forma una medida más 
confiable para monitorear el riesgo por inunda- 
ción en los 17 municipios del estado de Tabasco. 


Conclusiones 


El mapa de riesgo obtenido fue contrastado con 
los resultados históricos sobre desastres por 
inundación en los 17 municipios de Tabasco. 
Los resultados obtenidos fueron acordes con 
los registros de inundaciones en los años 1982, 
1989, 1995, 1998, 2002, 2003, y de 2005 a 2009. 
En el cuadro 10 se muestra el comparativo del 
índice de riesgo propuesto (IRP) y los registros 
de inundación en los años mencionados (IH). 

En algunos municipios no coincide, esto se 
puede deber a que no se encontraron los regis- 
tros de todos los municipios que presentaron 
problemas de inundación o que faltó considerar 
algunos otros factores que podrían aumentar el 
grado de vulnerabilidad o de costo. 

Finalmente, se puede concluir que en ge- 
neral los factores incluidos en la construcción 
del índice fueron correctos, así como el uso 
de componentes principales para obtener las 
ponderaciones adecuadas para cada uno de los 
factores considerados. 


Agradecimientos 


Los autores agradecen el apoyo de la LP-15 del Colegio 
de Postgraduados y de la LGAC-2 de la MII de la SEPI- 
UPIICSA-IPN para esta investigación, así como a la Comi- 
sión Nacional de Agua por la información sobre la precipi- 


tación de lluvias y los mapas de las regiones hidrográficas 


Referencias 


Arreguín-Cortés, F. I., Rubio-Gutiérrez, H., Domínguez- 
Mora, R., €: Luna-Cruz, E. (abril-junio de 2014). Análisis de 
las inundaciones en la planicie tabasqueña en el periodo 
1995-2010. Tecnología y Ciencias del Agua, 5(3), 5-32. 

Baró-Suárez, J. E., Díaz-Delgado, C., Calderón-Aragón, 
G., Esteller-Alberich, M., €££ Cadena-Vargas, E. (julio- 
septiembre de 2011). Costo más probable de daños por 
inundación en zonas habitacionales de México. Tecnología 
y Ciencias del Agua, 2(3), 201-218. 

Bayliss, A., £ Reed, D. (2001). The Use of Hisorical Data 
in Flood Frequency Estimation. Center for Ecology and 
Hidrology. Wallingford, England: MAFEF. 

Blaikie, P., Wisner, B., Cannon, T., €: Davis, L (1994). At Risk. 
Natural Hazards, People“S Vulnerability and Disasters. New 
York: Routledge. 

Campos-Aranda, D. F. (2001). Contraste de cinco métodos 
de ajuste de la distribución GVE en 31 registros históricos 
de eventos máximos anuales. Ingeniería Hidráulica en 
México, 7(2), 77-92. 

Cenapred (2006). Guía básica para la elaboración de atlas 
estatales y municipales de peligros y riesgos. México: Centro 
Nacional de Prevención de Desastres. 

Coles, S. (2001). An Introduction to Statistical Modeling of 
Extreme Values. London: Great Britain Springer-Verlag. 
Coles, S., £ Dixon, M. (1999). Likelihood-Based Inference 
for Extreme Values Models Extremes. Journal of the Royal 

Statistical Society, 2(1), 5-23. 

Coronel-Brizio, H. F., €r Llanos-Arias, J. (1996). Sobre el 
ajuste de la distribución Gamma a datos de precipitación 
pluvial. Ingeniería Hidráulica en México, 11(1), 15-20. 

Díaz-Delgado, C., Bá, K. M., € Trujillo-Flores, E. (1999). Las 
funciones Beta-Jacobi y Gamma-Laguerre como métodos 
de análisis de valores hidrológicos extremos. Caso de 
precipitaciones máximas anuales. Ingeniería Hidráulica en 
México, 14(2), 39-48. 

Duncan, W. (2002). Reinforcing Flood-Risk Estimation. 
Philosophical Transactions: Mathematical, Physical and 
Engineering Sciences, 360(1796), 1373-1387. DOI: 10.1098/ 
rsta.2002.1005. 

Fisher, R., € Tippett, L. (1928). Limiting Forms of the 


utilizadas. Frequency Distribution of the Largest or Smallest 
Cuadro 10. Comparativo del índice de riesgo propuesto contra los registros históricos. 

Mun. |Balancán | Cárdenas Centla Centro | Comalcalco | Cunduacán | E. Zapata | Huimanguillo | Jalapa 
IH Bajo Alto Muy Alto Alto Alto Alto Alto 
IRP Bajo Alto Alto Medio Alto Muy alto Alto 

Mun. Jalpa Macuspana | Nacajuca Paraíso Tacotalpa Tenosique 
IH Medio Medio Medio Medio Bajo Medio 
IRP Medio Medio Bajo Medio Medio Medio 


Tecnología y Ciencias del Agua, vol. VII, núm. 2, marzo-abril de 2016, pp. 99-114 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VII, núm. 2, marzo-abril de 2016, pp. 99-114 


Reyes-Olvera y Gutiérrez-González, Modelación del riesgo de inundaciones en el estado de Tabasco en el periodo 1961-2007 


Memeber of a Sample. Proceedings of the Cambridge 
Philosophical Society, 24, 180-190. 

Flowers-Cano, R. S., Jeffrey-Flowers, R., $: Rivera-Trejo, F. 
(septiembre-octubre de 2014). Evaluación de criterios 
de selección de modelos probabilísticos: validación con 
series de valores máximos simulados. Tecnología y Ciencias 
del Agua, 5(5), 189-197. 

Gnedenko, B. (1943). Sur la distribution limite du terme 
maximum d'une serie aleatoire. The Annals of Mathematics, 
44(3), 423-453. 

González-Camacho, J. M., Pérez-Rodríguez, P, €: Ruelle, 
P. (octubre-diciembre de 2011). Estimación de índices 
normalizados de lluvia mediante la distribución gamma 
generalizada extendida. Tecnología y Ciencias del Agua, 
2(4), 65-76. 

Green, C., Van Der Veen, A., Wiertra, E., €: Penning Rowsell, 
E. (1994). Vunerability Refined: Analysing Full Flood Impacts. 
London: Middlesex University Press. 

Gumbel, E. (1934). Les moments des distribution limite du 
terme maximum dúne serie aléatoire. Annals Mathematics, 
198(1), 141-143. 

Gumbel, E. (1958). Statistical of Extremes. New York: 
Columbia Univ. Press. 

Gutiérrez-González, E., Panteleeva, O. V., £ Córdoba-Lobo, 
V. M. (2012). Maximum Likelihood Estimators for the 
Parameters of a Distribution Extended Log-Loc-Scale 
Class. Far East Journal of Theoretical Statistics, 40(2), 127- 
143. 

Gutiérrez-González, E., Villaseñor-Alva, J. A., Panteleeva, 
O. V., £ Vaquera-Huerta, H. (2013). On Testing the Log- 
Gamma Distribution Hypothesis. Computational Statistics, 
28(6), 2761-2776. 

Herrera-Díaz, 1., Rodríguez-Cuevas, C., Couder-Castañeda, 
C., é€ Gasca-Tirado, J. R. (enero-febrero de 2015). 
Modelación numérica hidrodinámico-hidrológica en 
zonas de inundación con presencia de infraestructura. 
Tecnología y Ciencias del Agua, 6(1), 139-152. 

Hunt, J. (2002). Floods in a Changing Climate: A Review. 
Philosophical Transaction, 360(1796), 1531-1543. 

Jenkinson, A. (1955). The Frequency Distribution of the 
Annual Maximum (or Minimum) Values of Meteorological 
Elements. Quart. J. Roy. Meteo. Soc., 81, 158-171. 

Lehner, B., £ Dóll, P. (2001). Europe's Flood Today and in 
the Future. Model-Based Assessment of European Water 
Resources and Hidrology in the Face of Global Change. 
Center for Envioromental Sistems Research, University of 
Kassel, 5, 6.1-6.14. 

Merz, B., Thieken, A., € Gocht, M. (2007). Flood Risk 
Management in Europe: Innovation in Policy and Practice. S. 
S. Begum (Ed.). Berlin: Springer. 

Mileti, D. (1999). Disasters by Design. A Reassessment of 
Natural Hazard in the United States. Washington, DC: 
Joseph Henry Press. 


ISSN 0187-8336 


Ordaz, M. (1996). Algunos conceptos del análisis de riesgo. 
Cenapred, 14, 6-12. 

Rivera-Trejo, F., Soto-Cortés, G., € Barajas-Fernández, J. 
(octubre-diciembre de 2009). La inundación de 2007 en 
Tabasco, México: evolución de niveles de agua. Ingeniería 
Hidráulica en México, 24(4), 159-166. 

Rosbjerg, D., £ Madsen, H. (1996). The Role of Regional 
Information in Estimation of Extreme Point Rainfall. 
Atmos. Res., 42, 113-122. 

Smith, R. (2001). Flood Damage Estimation. A Review of 
Urban Stage-Damage Curves and Loss Function. Water 
S.A., 20, 231-238. 

Tapsell, S., Penning-Rowsell, S., €: Tunstall T. (2002). 
Vulnerability to Flooding: Health and Social Dimensions. 
Philosophical Transactions: Mathematical, Physical an 
Engineering Sciences, 360, 1511-1525. 

Tawatchi, T., £ Mohammed, F. (2005). Flood Hazard and 
Risk Analysis in the Southwest Region of Bangladesh. 
Hydrological Processes, 19, 2055-2069. 

Tawn, J. (1992). Estimating Probabilitis of Extreme Sea- 
Levels. Journal of the Royal Statistical Society, 41(1), 77-93. 
Uribe-Alcántara, E. M., Montes-León, M. A., € García-Celis, 
E. (abril-junio de 2010). Mapa Nacional de Índice de 

Inundación. Tecnología y Ciencias del Agua, 1(2), 73-85. 

Von Mises, R. (1936). La distribution de la plus grande de n 
valeurs. Revue Mathématique de lÚnion Interbalkanique, 1, 
141-160. 

Young, P. (2002). Advances in Real Time Flood Forecasting. 
Philosophical Transaction: Mathematical, Physical an 
Engineering Sciences, 360(1796), 1433-1450. 

Zucarelli, G. V. (noviembre-diciembre de 2013). Identificación 
de eventos hídricos extremos en la cuenca del río Paraná. 
Tecnología y Ciencias del Agua, 4(5), 181-187. 


Dirección institucional de los autores 


M.C. Ana Laura Reyes Olvera 


Colegio de Postgraduados 

Posgrado en Socioeconomía, Estadística e informática 
Campus Montecillo 

Km 36.5 carretera federal México-Texcoco 

56230 Montecillo, Texcoco, Estado de México, MÉxICO 
Teléfono: +52 (55) 5804 5900, extensión 1425 
mala_9860hotmail.com 


Dr. Eduardo Gutiérrez González 


Instituto Politécnico Nacional 

Sección de Estudios de Posgrado e Investigación 
Av. Te 950 Delegación Iztacalco Col. Granjas México 
08400, Ciudad de México, México 

Teléfono: +52 (55) 5624 2000, extensión 70276 
egutierrezgOipn.mx 


